Theoretical study of a localized quantum spin reversal by the sequential injection of 

spins in a spin quantum dot 
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This is a theoretical study of the reversal of a localized quantum spin induced by sequential 
injection of spins for a spin quantum dot that has a quantum spin. The system consists of "elec- 
trode/quantum well(QW)/dot/QW/electrode" junctions, in which the left QW has an energy level 
of conduction electrons with only up-spin. We consider a situation in which up-spin electrons are se- 
quentially injected from the left electrode into the dot through the QW and an exchange interaction 
acts between the electrons and the localized spin. To describe the sequentially injected electrons, 
we propose a simple method based on approximate solutions from the time-dependent Schrodinger 
equation. Using this method, it is shown that the spin reversal occurs when the right QW has energy 
levels of conduction electrons with only down-spin. In particular, the expression of the reversal time 
of a localized spin is derived and the upper and lower limits of the time are clearly expressed. This 
expression is expected to be useful for a rough estimation of the minimum relaxation time of the 
localized spin to achieve the reversal. We also obtain analytic expressions for the expectation value 
of the localized spin and the electrical current as a function of time. In addition, we found that 
a system with the non-magnetic right QW exhibits spin reversal or non-reversal depending on the 
exchange interaction. 

PACS numbers: 75.60. Jk,73. 63. Kv,85.75.-d 
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I. INTRODUCTION 

Magnetization reversal by spin injection (MRSI), in 
which the magnetization of a ferromagnet (FM) is re- 
versed by injecting a spin-polarized current into the 
FMjii^i^i^i^i^ili^i^ii^ is one of the most interesting topics in 
the field of spin-dependent transport and spin electronics 
for the following reasons: First, the reversal is induced by 
the current and not by the conventional method, such as 
the application of magnetic fields. Second, such a phe- 
nomenon is expected to have potential applications for 
the writing of data in the magnetic memory. 

Originally, the MRSI was theoretically predicted by 
Slonczewskii and Bergeri^ Their systems are FM/non- 
magnetic layer/FM junctions, in which the current fiows 
perpendicularly to the plane. Here, the magnetization of 
the FM is described by classical spins, and an exchange 
interaction acts between the magnetization and the elec- 
tron spiniii^ When the magnetization of one FM is free to 
change its direction and that of the other FM is pinned, 
the magnetization of the former can be changed to be- 
come either parallel or antiparallel to that of the latter 
depending on the direction of the current.^ In addition, 
some theoretical studies, in which the magnetization of 
the FM is described by the classical spins, have been re- 
ported since theni^i^ 

Quite recently, the MRSI has been experimentally 
observed in FM/non-magnetic layer/FM junctions^i^i^ 
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and pillars.— The experimental results have been often 
analyzed by using the theoretical model^ based on the 
classical spins. The magnetization of the FM appears to 
be fairly well described by the classical spins. 

On the other hand, we expect that magnetic materi- 
als will undergo transformation from classical spin sys- 
tems into quantum spin ones along with miniaturiza- 
tion towards high-density integration devices in the fu- 
ture. For example, such quantum spin systems are Mni2 
magnetic moleculesi ^^i^^'^'^i^'^'^^i^^ A Mni2 molecule pos- 
sesses an effective spin of S'=10 due to an antiferromag- 
netic interaction between the eight Mn'^+ {S=2) ions 
and the four Mn^+ (5=3/2) ions,— li^ii^ and it has a 
uniaxial anisotropy energy, — ji'lS'^, with D being an 
anisotropy constant with the magnitude of 0.7 K.— 
The anisotropy energy shows a bistable potential between 
Sz — ^Q and —10 states. Regarding a characteristic phe- 
nomenon due to this potential, the quantum tunneling of 
magnetization has been experimentally observed under 
the magnetic fields. This phenomenon has been ana- 
lyzed by using the quantum spin model with anisotropy 
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energy liiSi 

If the spin-polarized current can be injected into such 
quantum spin systems, the conventional magnetization 
reversal may be replaced by the localized quantum spin 
reversal. In the present condition, however, very few the- 
oretical studies for quantum spin reversal have been re- 
ported. Our primary focus is on what models bring about 
the quantum spin reversal and the length of the reversal 
time of the localized spin. The latter is important for 
the purpose of roughly estimating the minimum relax- 
ation time of the localized spin to achieve the reversal. 

In this paper, we examined the quantum spin 
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reversal induced by the sequential injection of 
spins for a spin quantum dot having the quantum 
spin. The system consists of "electrode/quantum 
well(QW)/dot/QW/electrode" junctions, where the left 
QW (L) has an energy level of conduction electrons 
with only up-spin. We considered a situation in which 
up-spin electrons were sequentially injected from the 
left electrode into the dot through the L. To describe 
the sequentially injected electrons, we first proposed a 
simple method based on the time-dependent Schrodinger 
equation. Using this method, we obtained an expression 
of the reversal time of the localized spin when the right 
QW (i?) had energy levels of conduction electrons with 
only down-spin. Furthermore, analytic expressions for 
the expectation value of the localized spin and the 
electrical current were obtained as a function of time. 
We also found the spin reversal occurred for the case 
of the specific exchange integral even when the R was 
non-magnetic. 

The present paper is organized as follows. In Sec. [TTl 
we present the model of the spin quantum dot and as- 
sumptions for the sequential injection of spins. In Sec. 
mil we provide a theoretical formulation; the wave func- 
tion, the electrical current, and the expectation value of 
localized spin are derived. In Sec. IIV A[ this theory is 
applied to the case of spin-polarized R, while Sec. lIVBl is 
the application to the case of non-magnetic R. In Sec. |Vl 
we make a proposal for a model with reversible switching 
and discuss the experimental aspects. Section IVII is the 
conclusion, and Sec. I VIII is the appendix, which includes 
information about the bias and gate voltages. 



and the R. The localized spin has a uniaxial anisotropy 
energy showing a bistable potential, — \D\S-^. The energy 
levels of the localized spin can be characterized by ■ It 
is also assumed that the magnetic easy axis of the local- 
ized spin is coUinear to that of the electron spins at the L 
and the R. In addition, the magnetic couplings between 
the localized spin and the QWs are negligibly small. 




5 = -S 5„=S 
Bistable potential 



FIG. 1: Schematic illustration of "electrode/L/spin quan- 
tum dot/_R/electrode" junctions. The L has an energy level 
of conduction electrons with up-spin, while the R has Nf en- 
ergy levels of conduction electrons with up-spin and Ni , with 
down-spin. The shaded area in the electrodes represents the 
region occupied by electrons. The bias voltage (the gate volt- 
age) is applied to the electrodes (the L and the R). In addi- 
tion, the equivalent circuit is shown in Fig. 1101 



II. MODEL 

A. Spin quantum dot 

In Fig. [U we show a system consisting of 
"electrode/L/spin quantum dot/i?/electrode" junctions, 
where the electrons flow from the left electrode to the 
right one under the bias voltage between the electrodes 
and furthermore the gate voltage is applied to the L and 
the R [see Sec. IVIIj . Here, the dot region behaves as a 
tunnel barrier, and tunnel barriers are set between the 
left electrode and the L and between the R and the right 
electrode. The L and the R form the QW; that is, the 
L has an energy level of conduction electrons with up- 
spin and becomes a spin filter to inject only the up-spin 
electrons, while the R has N^i energy levels of conduction 
electrons with up-spin and N^^ with down-spin. Those 
levels in the R are introduced to accept electrons ex- 
hibiting the elastic and the inelastic transports [see Sec. 
HTBl iv)]. The left electrode (right electrode) IS a non- 
magnet or ferromagnet, in which the direction of the spin 
polarization is the same as that of the L (R). The shaded 
area in the electrodes represents the region occupied by 
electrons. In the dot, the quantum spin S={Sx, Sy, Sz) 
with S' > 1 is localized and is weakly coupled to the L 



B. Assumptions 

On condition that the initial state of the localized spin 
is Sz=—S, up-spin electrons are sequentially injected 
from the left electrode into the dot.-l^ To describe the 
sequentially injected electrons, we have the following as- 
sumptionSf^ 

(i) This system exhibits single electron tunneling 
(SET),^^ in which the up-spin electrons are injected 
from the left electrode into the dot one by one un- 
der specific bias and gate voltages [see Sec. IVII| . 
Concretely speaking, when the probability density 
of the electron becomes at the L and 1 at the R, 
the electron moves into the right electrode and the 
subsequent electron is injected into the L. 

The situation in which the mth electron runs in 
the dot is named as the mth process. When the 
time period of the mth process is represented by 
At(™\ the finish time of the mth process, tm, is 
defined by t„ = J2Z'=i At^"'). Here, to=0 is the 
initial time. 

(ii) When the transport of the electron is caused by 
interaction V between the L and the R, V acts 
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while each electron runs in the L-dot-R system. On 
the other hand, V is switched off for the electron 
being outside of this system. The interaction for 
the mth injected electron at time T, V^™'>{T), is 
written by, 



1/(™)(T) = 



V, for t„,-i <T< t,n, 
0, for T < t,n-i, tm < T, 



(1) 



[see Fig. H]. As found from this y(™)(T), each 
process is independent. 

The interaction V is given by a transmission 
term including the exchange interaction between 
the electron and the localized spinji^i22iiSli2^ This 
term is obtained within the second-order perturba- 
tion theory based on the weak couplings between 
the dot and the QWs. The L, R, and dot parts 
are described by an unperturbed Hamiltonian, in 
which the on-site Coulomb energy is considered in 
the dot. The couplings between the dot and the 
QWs are perturbation terms. The resultant V is 
now expanded to a more general system, namely, 
the anisotropic exchange interactionf22. The expres- 
sion is written by. 



V = Vo 



nT = l 

nl = l 

N^ 

nT = l 



(2) 



Here, Cj (c^) is the annihilation (creation) opera- 
tor of an electron of j=L or Rnai where the suffix 
L (Rna) denotes the level of conduction electrons 
with up-spin of the L (the nth level of conduction 
electrons with a a spin of the R). Furthermore, Vq 
denotes the coefficient for direct tunneling, while 
J± (Jz) is that for tunneling with a transverse (lon- 
gitudinal) exchange interaction between the elec- 
tron and the localized spin. From now on, we call 
J± (Jz) the exchange integral. 

(iii) The localized spin canted by the electron is not 
rapidly relaxed and interacts with the subsequently 
injected electron. 

(iv) The system has the relation of |£'l — Eji^J <$:.\V\, 
where Ej^ {Er„„) denotes the energy level of con- 
duction electrons with the up-spin of the L (the 
nth energy level of conduction electrons with the 
a spin of the R). The energy levels contain the 
potential due to the gate voltage. Here, the cen- 
tral region in En^^ is located in the vicinity of E^. 
The band width inEa^^ is about 2\D\{2S-1), with 
\D\ [S^ -{S- 1)2] = |L>|(25 - 1), where \D\{2S ~ 1) 



corresponds to the maximum energy that the con- 
duction electron absorbs from the localized spin or 
gives to the spin for the case of the energy conser- 
vation. Consequently, \El — Eji^^\ <^ \V\ corre- 
sponds to \D\ <^ \V\. Note that the assumption 
greatly simplifies the calculation [see Sec. IHI A| . 

As mentioned in Sec. IHI Al [also Ref. 24], the 
above assumption leads to the relation of At^™^ <C 
h/\EL — Eii^^^\, where h is the Planck constant h 
divided by 2-k. On the basis of the relation, we 
investigate the motion of the electron within the 
L-doi-R system using the wave function obtained 
from the time-dependent Schrodingcr equation. In 
contrast, the Fermi's golden rule which is a method 
to investigate the transport property, is applicable 
for the case of At'™) ::^h/\EL- Er^^\. 

We also note that in a short time period such as 
At*"), the total energy becomes uncertain accord- 
ing to the uncertainty relation between time and 
energy, although the energy can be certainly con- 
served after enough time has passed. The electrons 
then exhibit the elastic and the inelastic transports. 
Such electrons are accepted by some energy levels 
in the R. 



V 



At 



(m) 



t.,-. t„ 



FIG. 2: Interaction for the mth injected electron in the mth 
process, V^"^\T). Each process is independent. 



III. THEORETICAL FORMULATION 

A. Wave function 

In order to obtain the wave function of the mth pro- 
cess, |$(™)(i)) with t = T-tm-i and < i < Ai("), we 
solve the time-dependent Schrodingcr equation. 



^ 7^(")(f)|$(")(f)), 



with 



(m)/ 



(3) 



(4) 



where Tip is the Hamiltonian for the electronic states of 
L, i?, and the localized spin state. The wave function 
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for < t < At^") is expressed as, 



|$M(t)) 



s 

J s, = -s 



(m) 



We- 



with 



(5) 



(6) 



forj=Lor i?„^, a.TLdLEL.,s,^EL-\D\Sl, Er^^^s,=Er^^- 
jDjS'^. The state \L {Rna),Sz) means that the electron 
exists at the level of the L (at the nth level of the a spin 
of the R) and the localized spin has Sz- The equation 
to determine the coefficient a^™'' (t) for < t < Ai^™) is 
given by, 



fe si = -s 



0-,^.|y|fc,^;)45,(t) 



xe 



We obtain equations for 



(m) 



and a;^^";^5^+i(t) using + l\V\L,Sz ) = 

(L, ^,|y|i?„x, Sz + 1) = J±y/{S - Sz){S + Sz + 1) = Fs^ 
and {R„^,Sz\V\L,Sz) = {L, Sz\V\R^^, Sz) = Vo+JzSz = 
Gs^ ■ Their equations are summarized as follows: 



it) 



,(m) 



F: 



(8) 



^':!,5.W = -^^ail(^)e-^«"T>-.*/^ (10) 



G 



with 



ft\t) 



9s 



G, 



AT, 



ni = l 



fi2 



nT = l 



A_R„^,s^ — El^s^ — Er^^^sz+Ii and A^^^^s^ = E^.s^ — 
Er^^.Sz, where the sum Y^ni (Snt) excludes ni (nj) 
giving Afl„^,s^=0 (A/j„^,s^=0). 

According to \El - Er^J < |Jj.|, \ Jz\, and |Vb| of 
Sec. Ill Br iv). the present system has the relation of 
\Ar^„sJFsA « 1 and \Ar^^^sJGsA < 1^ The time 

average of \f^\t) \ and |55"''(i)| then becomes negligibly 
small compared to that of the magnitude of the first and 
second terms on the right-hand side of Eq. ([5]) because 
the sum in the complex plane is done in addition to this 



relation and |a^| g +i(i)| and |a^"!j g {t)\ are 1 at max- 
imum. We hence obtain (t) from an approximate 
expression by omitting fg^\t) and gs"\t) from Eq. ([S]). 

The assumption for a^™J (t) of Sec. IIIBT i) is concretely 
defined as follows: under the voltages described in Sees. 
IIIBf i) and I VIII as soon as the local probability density at 
the L of the mth injected electron becomes zero, the mth 
electron goes out of the R, and the potential for the mth 
electron is then switched off: (T) = for T > t^. 

Next, the subsequent (m -|- l)th electron moves from the 
left electrode to the L. The initial condition of the mth 
process is therefore given by, 



Am) 



(m) 



.(0) = 



J2 \{Rn^,Sz\<f'^"'-'HAt("'-'^)) 



nT = l 



+ J2 |(i?ni,5,|cI>('"-l)(Ai("-l))) 

nl = l 
Am) 



41(0) 



Here, &<j'"'' is the initial amplitude of the mth process with 

< b^™'' < 1. It is represented by using the probability 
amplitude specified by the final state of the localized spin 
of the (m — l)th process according to Sec. IIIBr iii). In 

particular, b^;}g—l is apparent due to Sz=—S at T=0. 
Furthermore, d^"2 (0)=0 originates from V^"^\T)=0 for 
T<tm-i of Eq.^dl]). 

As a result, the coefficients are obtained as follows: 



= 0. 



1/2 

(12) 
(13) 



A":iit) 



6^"^ cosies J), 



Am) 



Am) 



iFs^h 



(m) 



(m) 



■sin(ris,i), 



(14) 
(15) 

(16) 



with 



(17) 



(m) 



for < < < At(™). Here, a 

are obtained by using a^l(i) and Eqs. (O and (fTU)) . 
where AR^^,s^t/h and AR^^^s^t/h are taken as because 
t is considered to be, at most, several times larger than 
l/^ls^ in this study^ and further \Ar^^^sJFs^ \ < 1 and 
\Ar^^^S^/GsJ <C 1 are introduced as stated earher. 

In fact, the coefficients of Eqs. ((HI), and (fTB)) 

correspond to exact expressions for a system of D=0 



and Ai^„^,s^=A/ 



=0 [El 



Ep 



Er^ ) obtained 
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under the condition that the initial state of the local- 
ized spin is set to be Sz=—S. The use of the coefB- 
cients, however, is considered to be valid for the study 
of the qualitative properties of the present system with 
\El-Er^J < |Ji|, |J,|, and \Vo\. 

We substitute Eqs. and HH) into Eq. 

The wave function |<I>(™)(t)) for < t < At'™) is finally 
written as, 

S 

s^=-s 

S 

Ni S-1 

(18) 



(m) _ "iTiW Am) 



With AHft) 



, (f) - 



a 



(m) 



At) 



normalization factor, 



, where a^"^\t) is the 



E haW| E 



1/2 



(19) 



The local probability densities at the L and the R are 
given by. 



s^=-s 



(20) 



E KIa+iw', (21) 



respectively, with £ 



J2s, = -S 1-^' "^z) (-^' '^■s 



and n = J2n{=lJ2s,=-s\RnhSz){Rn],S;, 



B. Number of electrons in the R and the right 
electrode 



On the basis of Eq. pT|) . we obtain the number of 
electrons in the R and the right electrode for the mth 



process, ni"^\t). The expression is written by, 

s 



s^=-s 



(22) 



for < i < Ai('"). The first term in the right- 
hand side represents the number of electrons which have 
moved to the right electrode. The second and third 
terms are ($(™) (^)|7^|$('") (0) ofEq. Furthermore, 
n^r^\t)-{m^l) ($(™)(^)|7^|$(™)(^))] is regarded as 
the expectation value of the position of the mth injected 
electron when the positions of the L and the R are set to 
be and 1, respectively. 



C. Electrical current 

The electrical current of the mth electron, I^"^^{t), is 
defined by /(™) (t)=-eneS'dw(™HO> where w(™)(<) is the 
velocity of the mth electron, e (>0) is the electric charge, 
rig is the number density of electrons, and Sd is the cross- 
sectional area of the dot. The velocity, u(")(t), is ob- 
tained from the time differential of the expectation value 
of the position of the mth injected electron. When the 
positions of the L and the R are defined as and I, re- 
spectively, /(™)(f) is written as follows: 



/M(t) = _ene5d/^(<i>("^^)|7^|<i'<"^(^)> 
at 



S,l 



(m) 



(23) 



for < t < Ai(™). 



D. Expectation value of the localized spin 

The expectation value of the localized spin of the mth 
process, {Si"'\t)), is obtained by using Eq. ([15]). The 
expression is given by, 

(5(™)(<)) = ($(™)(t)|5,|$('")(0) 



s^=-s s^=-s 



s,=~s 
for < i < At(™). 



(24) 
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IV. APPLICATIONS 
A. Spin-polarized R with only down-spin levels 

The spin reversal occurs when the R has energy lev- 
els of conduction electrons with only down-spin. In this 
system, V is described by only the second term in Eq. 



because of Ni=0. Furthermore, Gs^, sA^^' 



A 



(m) 



(t) are set to be 0. 



Q 



>>0.5 



o 




0.25 

FIG. 3; Time t dependence of the partial probabil- 
ity density of the mth process. Here, \{L, Sz\$'-"^\t))\^ 

iT,ni=i\iRni,S, + l\^^"'\t))f) is shown by a solid (dashed) 
curve. 



We obtain for < i < At^'") using Eqs. dH]), 

(Uni), CH]), and 6^"^= (5„,s^+5+i which is Eq. ^ for this 
system. The expression is written by, 

= cos(a;s.i)e-*^^'^^*/'"lL,^,) 
iJ± 



nl = l 



3in(^S^t)e-'^«"i'^^+^*/''|i?„X,5. + l), 

(25) 



with 



\J±\^NaS-S,){S + S, + l) 



(26) 



and m=S^ 5" -I- 1 for S^=-S -5-1. From Eq 
the partial probability density is obtained as follows: 

\{L,S,\^^"'\t))\^ ^cos\u;sJ), (27) 

I 2 

^ |(i?„X,^. + l|*^"^W) =sin2(c.s.t), (28) 



with m=Sz 



5 + 1 for S^=-S - S - 1. For m = 



5"^ -t- S* -I- 1, \{L,S^\¥"'\t))\ decreases from 1 to 

and J2ni=i \ {Rni,S^ + l|$(")(i))|^ increases from to 
1 with increasing t [see Fig. [3]. This behavior means 



that the electron transfers from the L to the R with the 
spin flip from up to down and the localized spin changes 
from Sz to Sz + I- Since the mth process is finished at 
t = tt/{2ujsJ, Ai(™) is given by. 



2^s. i\J±\^Ni{S^ Sz){S + Sz + l) 



,(29) 



with m = Sz + S + I for Sz——S - S — 1. Here, m = 
Sz + S +1 for Sz^-S - S* - 1 represents m=l, 2, 3, 4, 
• • •, 2S. The spin reversal is thereby achieved after the 
(25)th process, where 25 also corresponds to the total 
number of injected electrons. 

Using Eq. ([^5]) , we can obtain the reversal time of the 
localized spin, tn, which is the total time to transform 
Sz from -S to S. When At(") with m = Sz + S + 1 of 
Eq. (|29p is renamed as Atg^ , tji is simply written by, 



S-l 

= y = 

s~is 4 l^^l V^liS- Sz){S + Sz + l) 



h 



(30) 



Figure m shows the S dependence of tii\J±\y/Ni/h. At 
S=l,tii takes ft.V2/(4| J_i|y]V|). On the condition that 
\J±\y^ Ni/h is a constant, tj^ monotonically increases 
with increasing S. This behavior is explained by con- 
sidering that the spin reversal is caused by injecting 25* 
up-spin electrons into the dot. In the classical spin limit 
of S ^ oo, tji comes close to /i7r/(4| J^|-y]Vj") which is 

obtained by replacing J2s =-s with Muis^oo J^g^ dSz- 
We thus find that the tn of systems with S > 1 has the 
following relation. 



hV2 



4|J, 



<tB. < 



hn 



(31) 



Furthermore, since this is derived under Sec. 
IIIBf iii). we emphasize that the spin reversal is re- 
alized when the relaxation time of the localized 
spin, r, satisfies the relation of r > tj^ — 

Y.l:Lsh/ [^\J±WNi{s-Sz){s + Sz + i) . 

As an application, we consider a system with S=h, 
iV|=50, and |J_l|=0.001 eV. The reversal time, Ir, is 
evaluated to be 3.3 x 10" s by using Eq. The re- 

lation of T > 3.3 X 10""'^^ s is essential for an experimental 
observation of this spin reversal. 

The upper panel of Fig. [5] shows the T dependence 
of nr^\T — tm-i) with TO=1 - 10. The number of elec- 
trons in the R and the right electrode, nf^\T — tm-i), 
increases from to 10 with increasing T and becomes 10 
aX T = tn = 3.3 X 10~^^ s. The reversal is eventually 
achieved by injecting 10 (=25) electrons into the dot. 
Furthermore, it is noted that the slope becomes zero at 
tm with m—\ - 9. 
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0.4 ° 



0.2 



° : tR\JJN\'^/h 
---: 71/4 



10° lO' 10^ 10^ 10"^ 
S 

FIG. 4: Spin quantum number S dependence of the reversal 
time of the localized spin, Ir. Here, tR\J±\y/W[/h takes \/2/4 
at S=l, while it comes close to 7r/4 in the classical spin limit 
of 5 ^ oo. 



The middle panel of Fig. [5] shows the T dependence 
of /(™)(r - i„_i)/(-ene5dO- The quantity /(™)(r - 
tm-i) / {—eneSdl) exhibits a single peak for each process. 
The peak takes the largest value for the 5th and 6th pro- 
cesses, where the respective initial states have ^2 = — 1 
and 0. This behavior reflects the fact that ujs^ achieves 
its maximum at 5^ = — 1 and for S=5 [see also Eq. 

mi 

The lower panel of Fig. [5] shows the T dependence 
of {Si"^\T — tm-i))- The expectation value of the lo- 
calized spin {Si"^\T — tjn-i)) increases from —5 to 5 
with increasing T following the same pattern as that 
of ni"^\T — tm-i)- It finally becomes 5 at T = t^j = 
3.3 X 10^^^ s; that is, the spin reversal is completed at 
this tfj. 



B. Non-magnetic R 

In the case of non-magnetic R with N^=Ni, this sys- 
tem exhibits spin reversal or non-reversal depending on 
the exchange integral. 

Regarding the spin reversal, all injected electrons can 
transfer from the L to the R; that is, the local probability 
density at the L becomes zero at a certain time for each 
process [see Sec. IllBr i)]. We therefore pay attention to 
the local probability density at the L of the mth process, 



s 

E 



cos^insj), (32) 



which is obtained from Eqs. ([H]) and ([^0]) . Equation 
([32)1 becomes zero at t = At^™), when At^™) satisfies the 
following relation. 
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FIG. 5: Time T dependence of the number of electrons in the 
R and the right electrode, nl"^\T — tm-i)', electrical current, 
_ and expectation value of the localized spin, 

{S^^\t - tm-i)); m=l - 10 for the system of S=5, N^^bO, 
and |Jx|=0.001 eV. Upper panel: nl""' (T - Middle 
panel: I^'^'> {T ~ tm-i)/{~eneSJ). Lower panel: {Si'^^T- 

tm-l))- 



for any giving finite 65^ 
counting number. 



Here, nl"' (> 0) is the 



In the first process with 6^5=1, Eq. 
cos^(ri_si), and it certainly takes zero at t 
7r/(2ri_s)- In the mth process with m > 2, 6^'"'' is finite 
for any Sz ■ Whether Eq. ((33|) is satisfied for m > 2 or not 
depends on the parameter set of the exchange integral. 

We first investigate the parameter set satisfying Eq. 
P3p in order to find the spin reversal. For a system with 
S—2 and iV|=iV|=20, the parameter set is, for example, 
^0=0.001 cV, J^/V^)=12/31, JJVq = 2711/31. In Table 
m we summarize b'^p\ n"g^\ and At^™) for each m. The 

first process with b''ll=l has At^^) = 7r/(2f^_2) = 4.78 x 
10~^^ s, and the mth process with m > 2 has A<(™) — 
37r/(2n_2) = 1.43 X lO'^^ s. 

The upper panel of Fig. [6]and the upper panel of Fig. [7| 
show the T dependence of — t,„_i) with the above 

parameter set in the regions ofO<T<8xlO~^^s and 
< T < 1 X IQ-" s, respectively. Here, T = 8 x 10"" 
s (T = 1 X 10~^^ s) corresponds to the time in the 
57th (8th) process. The number of electrons in the R 



becomes 
At(i) = 

(™) V 



(33) 
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FIG. 6: Time T dependence of n^' (T-tm-i) and {Sr'' (T- 
tm-i)) with m=l - 57 for the system of S=2 and N^=Ni=20. 
The parameter set of 1/ is l/o=0.001 eV, Jx/Vo= 12/31, and 
Jz/Vo = 2\/TT/31. Upper panel; ni"^\T — tm-i)- Lower 
panel: {S*™'(r - 



and the right electrode, rir {T — tm-i), slightly oscil- 
lates in each process and clearly increases with increas- 
ing T. In the middle panel of Fig. [71 it is shown that 
ji'^) (^T—tjn-i)/ {—eUeSdl) oscillates between positive and 
negative values 

The lower panel of Fig. [Hland the lower panel of Fig. [7] 
show the T dependence of (S'1™-'(T— tm_i)) in the respec- 
tive regions of T. The expectation value of the localized 
spin, {si"^\T — tm-i)), has a slight oscillation in each 
process as well. The spin reversal is almost realized in 
the vicinity ofr = 8.0xl0^^^s, and it is also confirmed 
from h^f^ = 0.99 in Table HI 

As the spin non-reversal case for 5=2 and N^^~N ^=20, 
we choose, for example, the parameter set of Vb=0.001 
eV, J_i/Vb=l, Jz/Vq=2, which does not have t to satisfy 
Eq. (1221) for Sz = -2,-1,0,1,2. In the first process 

with 6L^^=1, Eq. dMl) becomes \{^^''\t)\C\^'^^Ht))\^ oc 
cos^(rj_2t) and it takes zero at < = 7r/(2f7^2)=6.4x lO"^** 
s. In fact, the first injected electron is transferable as seen 
from the upper panel of Fig. [H In the second process, 
Eq. (122) always has finite values. The second injected 
electron is trapped in the dot and it oscillates between the 
L and the R as shown in the upper and middle panels of 
Fig. [21 The localized spin also oscillates between 5*2=- 2 
and 5*2=- 1 [see the lower panel of Fig. [5]. 

In addition, we report that qualitative behaviors iden- 
tical to those noted above are found for some systems 
(some S"s). 
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FIG. 7: Time T dependence of n['"'(r - /(™)(r - 

tm-i) / {—erieSdl) , and {Si"^\T — tm-i)) with m=l - 8 for the 
system of 5=2 and Ni—N^=20. The parameter set of V is the 

same as that of Fig. 11 Upper panel: n,'.'"'(r-f„_i). Middle 

panel: I^"'\T - trn-i)/i-en,SJ). Lower panel: {Si"'\T - 

tm~i)). Here, n,'"' (T - and {Si^'^T - tm-i)) are the 

respective ones on the larger scale of Fig. [6l 



V. COMMENTS 

What follows is a comment about an application of 
the system with the non-magnetic R. According to the 
present study, the system exhibits a spin reversal or non- 
reversal depending on the exchange integral. The ex- 
change integral is originally described by the orbital en- 
ergy level and the on-site Coulomb energy in the dot 
and transfer integrals between the dot and the QWs^ 
Thereby, by controlling the orbital energy level with 
the application of the gate voltage, this system may be 
switchable between the spin reversal and the spin non- 
reversal. 

Secondly, we propose a model with reversible switching 
between Sz=—S and S [see Fig. [3]. The L has energy 
levels of conduction electrons with up-spin, in which an 
energy level is located at Ei^l and the lowest of the iV| 
energy levels is found at £'2,l ■ The R has levels of conduc- 
tion electrons with down-spin, in which a level is located 
at £'2,fl and the lowest of the levels lies at Ei^n. Here, 
tunneling probabilities between the level at Ei^^ of the L 
and that at i?2,fl of the R and between the levels in the 
vicinity of E2.L of the L and those in the vicinity of Ei^n 
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fog™'; count- 



TABLE I: Initial amplitude of the mth process, 
ing number of the mth process, rig'"'; and time period of 
the mth process. At'""' for the system of 5=2, A/'|=iV|=20, 
^0=0.001 eV, Jx/Vb=12/31, and J^/I/q = 2^/31. 
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of the R are assumed to be negligibly small. The energy 
of the highest occupied state (i.e. the Fermi level) of the 
right electrode is changed by applying the bias voltage, 
while that of the left electrode is fixed between Ei^l and 
E2,L- It should be noted that some electrons lying be- 
tween the energy of the highest occupied state of the left 
electrode and that of the right electrode contribute to the 
transport. 

We now consider the SET regime as described in Sec. 
IVIII When the bias voltage is applied to the right elec- 
trode as shown in Fig. [5] (a), the up-spin electrons are 
injected from the L into the dot and the reversal from 
Sz=—S to S can be induced through the interaction of 
Eq. ^ with N-f—O. On the other hand, in the vicin- 
ity of £'2,_R, the interaction between the electron and the 
localized spin is written by. 



V = J\ 



E 

«T=i 



.^4 



(34) 



where the suffix (L„| ) denotes the level of conduction 
electrons with the down-spin of the R (the nth level of 
conduction electrons with the up-spin of the L), and 
is the transverse exchange integral between the localized 
spin and the electron in this energy region. When the 
bias voltage is applied to the right electrode as shown 
in Fig. [HI (b), the down-spin electrons are injected from 



CO 



I 




8; Time T dependence of n'T^T - t™-i), j('")(r 



FIG. 

tm-i)/{-en^SJ), and {sITHt - t,n^i)) with m=l, 2 for the 
system of S—2 and Ni—Ni=20. The parameter set of V 
is Vo^O.OOl eV, J±/Vo=l, and J^/Vo = 2. Upper panel: 
' ni"^\T -tm-i). Middle panel: (T - t,„_i)/(-eneS'd/). 
Lower panel: {Si"^\T — tm-i))- 



the R into the dot and the reversal from Sz=S to —S is 
achieved through the interaction of Eq. (p4| . 

Finally, we discuss the experimental aspects from the 
viewpoint of the magnitude of current density. In the 
case of S—5, 7V|=50, and |Jj^|=0.001 cV, the required 
maximum current density for the spin reversal, j,nax, 
is evaluated to be j,„^^= I^^^/Sd ^ 4.2 x 10® A/cm^, 
where the maximum current J^^x is given by I^^^ = 
(5.9x 10^3) X (1.6x10-19) ~ 9.4 x 10"^ A as shown in the 
middle panel of Fig. [SJ with neSdl=l and e~1.6 x 10^^^ 
C. Furthermore, Sd is set to be S'd=(1.5 x 10^^)^ cm^, 
which is obtained by assuming that the molecule with 
iS'=5 has a cubic structure that is about 1.5 nm per side. 
This length is roughly estimated using examples from a 
similar molecule, Mni2 with S=10}^ In terms of sustain- 
ability, the abovementioned j^^^^ is high for typical ferro- 
magnetic tunnel junctions, such as CoFeB/MgO/CoFeB 
junctions with critical current densities for the MRSI of 
7.8 X 10^ A/cm^ and 1.2 x lO'^ A/cnP.'^'^ On the other 
hand, it has been recently reported that some carbon 
nanotubes can individually carry currents with density 
exceeding 10^ A/cm^ i^^i^^ We therefore anticipate that 
appropriate use of carbon nanotubes encapsulating mag- 
netic molecules'^" will lead to the achievement of the spin 
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(b) 

FIG. 9; A model with reversible switching between Sz~~S 
and 5*. (a) The reversal from Sz=~S to 5* is induced by 
injecting the up-spin electrons from the L into the dot. (b) 
The reversal from Sz~S to —S is achieved by injecting the 
down-spin electrons from the R into the dot. 



reversal. In addition, we strongly believe that the current 
density can be tuned by controlling the exchange integral 
and iVf (or Ni), as found from Eq. (^5)1 . Namely, the 
current density decreases with reducing the magnitude of 
the exchange integral and N-^ (or 7V|). 

VI. CONCLUSION 

We studied the localized quantum spin reversal for a 
case in which up-spin electrons were sequentially injected 
into the spin quantum dot. To describe the sequentially 
injected electrons, we first made assumptions about the 
sequential injection of spins and then proposed a simple 
method based on approximate solutions from the time- 
dependent Schrodinger equation. This method was ap- 
plied to systems with the following R; that is, (a) the 
spin-polarized R having energy levels of conduction elec- 
trons with only down-spin, and (b) the non-magnetic R. 
The respective results are summarized below. 

(a) The system exhibited a spin reversal. In particular. 



we derived the expression of the reversal time of the 
localized spin and explicitly showed the upper and 
lower limits of this time for systems with > 1 . In 
addition, analytic expressions for the expectation 
value of the localized spin and the electrical current 
were obtained as a function of time. 

(b) The system of 5=2 and N^=N i=20 exhibited the 
spin reversal or non-reversal depending on the ex- 
change integral. In the spin non-reversal case, the 
second injected electron was trapped in the dot, al- 
though all injected electrons transferred from the 
L to the R for the spin reversal case. 

We expect that the above-described phenomena will 
be observed with the advancement of experimental tech- 
niques and will be utilized in spin electronics devices in 
the future. 
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VII. APPENDIX: BIAS AND GATE VOLTAGES 
FOR THE SET 

The SET region in the plane of the bias and gate volt- 
ages for the model of Fig. [T] is found. We investigate 
the energy change AE which is obtained by subtracting 
the energy of the initial charge state from that of the 
final charge state. The change of state really occurs for 
AE < 0. 
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FIG. 10; Equivalent circuit for "electrode/L/spin quantum 
dot/J?/electrode" junctions. The tunnel barrier between the 
left electrode and the L (between the R and the right elec- 
trode) is named 1 (2). Each symbol composed of two rectan- 
gles represents both the condenser and the tunnel junctions 
which enable electron tunneling to take place. Each condenser 
of the gate is assumed to have a thick barrier that suppresses 
the electron tunneling. 

The equivalent circuit for this model is shown in Fig. 
[TUl The bias voltage of the left electrode (the right elec- 
trode) is represented by vl (vr), and the voltage of the 
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gate electrode is Vg. The capacitances of junctions of 1, 
2, the dot, and the gate are Ci, C2, C^, and Cg, respec- 
tively. The electric charge of the respective junctions is 
(i=l-4, d). 

When the initial charge at the L (R) is represented by 
QL {qn), the final charge at the L (R) is (9^)1 and 
the change of the charge number in the left electrode 
(right electrode) is An^e (An^e), the energy change 
^^9i,<ZH:9i:9kAni,e,A«^., is given by, 

(35) 

where AU denotes the change of the electrostatic energy, 
and Wi {W2) represents the work done by the voltage of 
the left electrode (the right electrode) and W3 (W4) is 
that of the left gate electrode (the right gate electrode). 
These expressions are written as follows: 



AU ^U' - [/, 



^i^L-Vgf + ^( 



2 



2 

a 

2 ''^'^ ' 2 

Wi = -AqiVL - eAuLeVL 

W2 = Aq2Vji + eAuReVR, 
W3 = -AqsVg, 
W4 = ~Aq4Vg, 
with 



(36) 
(37) 



^gf + ^{4>'R~Vg)\ 



(38) 

(39) 
(40) 
(41) 
(42) 



Agi 


= Ci(^ 




VL)-Ci{cj)L- 


vl), 


(43) 


Aq2 


= C2[vr - ( 


1>'r)-C2{vr~ 


<1>r), 


(44) 


Aq3 


= Cg{c 




^g) - Cg{(t)L - 




(45) 


Aqi 


-Cg{c 








(46) 




= aqR - 




+ 7, 




(47) 


4'R = 


= CqR^ 


- aqL 


+ V, 




(48) 


^y- 


= (^<i'r - 




+ 7, 




(49) 




-- Cq'R ^ 




+ '7, 




(50) 




Cd 








(51) 


a = 










P = 


— (C2 




+ Cg), 




(52) 



X 



[Cd{C2VR + CgVg) 



(53) 
(54) 



HC2+Cd + Cg){ClVL+CgVg)] 

(:^^{C4 + Ca + Cg), 

?7 = Y \Cd{ClVL+CgVg) 

+ {Cl+Cd + Cg){C2VR + CgVg)], (55) 

X^{Ci+Cd + Cg){C2 + Cd + Cg) - Cl (56) 



where e (>0) is the electric charge. The works, 
—eAriL^VL of Eq. ([5^ and eAuR^VR of Eq. (^0]) rep- 
resent that the electrodes are provided with the changed 
portion of the charge due to the tunnel electron. Fur- 
ther, (pL {4'r) is the electrostatic potential of the L {R), 
and Eqs. ([T7)) and (|48p are obtained from the following 
equations. 



qR 



-qd 

-92 



qi + qs, 
qd + qi, 



(57) 
(58) 



with qi^Ci{(j)L-VL), q2=C2{vR- (j)R), q3^Cg{(t>L - Vg), 
qi^Cg{cj)R - Vg), qd^Cd{(j)R - (f>L)- The bias volt- 
age vl {vr) is now set to be ~v/2 {v/2). When 



AE, 



=0, the relation between v and 



Vg is given by. 



- [\[di - ds) + iy{d2 + ^3)] Vg - dj 
K{di - ds) + ^{d2 + da) - ^iAuRe + AnLe)e 



,(59) 



with 



di = (Cl + Cg) [a{q'R - qR) + Piq'L ~ qL)] , (60) 
d2 = {C2 + Cg) [({q'R - qR) + aiq'^ - q^)] , (61) 
d3 = Cd [{( - a){q'R - qR) + (a - /3)((7i - qL)] , 

(62) 

C2+C„ 



C?4 



2 

Ci + C, 
2 

Cd 



[^(CqR + aqL)^ + {Cq'R + aq'L)^ 



[iaq'jj + Pq'^)^ ~ {aqR + PqL 



^{[(C-«)9k + («-/3)9L]' 
[iC-a)qR + ia-p)qL]^}, 



1^ — ■^j^{C2Cd — C1C2 — CiCd — CiCg), 



1 

X 

1 

2X 
1 



X=^iC2Cg + 2CdCg+C^g), 



(C1C2 -I- C2Cd + C2Cg — CiCd), 



X 



{C^Cg + 2CdCg+Cl). 



(63) 
(64) 
(65) 
(66) 
(67) 



This model exhibits the SET under the specific v and 
Vg, which satisfy A£;o,o,-e,o4,o < 0, Ai;_e,o,o -e,o,o < 
0, A£'o_e.o, 0,0,1 < 0, Ai;o.-e -e.-e,i.o > 0, and 
A£;o,o,o,+e,o,i > 0. In Fig. [TTl the SET region for a 
system with Ci=C2, Cc;=10Ci, and Cg=0.1Ci is shown 
by the shaded area. 

We furthermore describe the present transport prop- 
erty; that is, when the probability density of the injected 
electron becomes at the L and 1 at the R, the elec- 
tron moves into the right electrode, while the electron 
cannot go to the right electrode if the probability den- 
sity is finite for both the L and the R. This property is 
based on the following two assumptions: First, there is 
no direct transfer integral between the L and the right 
electrode. Second, the wave function extended over both 
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the R and the right electrode can be ignored, although 
that extended over both the L and the R is taken into 
account, under the condition that the coupling between 
the R and the right electrode is much smaller than that 
between the L and the R, i.e. Eq. 
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FIG. 11: SET region in the bias (u) and gate (vg) volt- 
age plane for the equivalent circuit with vl=—v/2, vr=v/2, 
Ci=C2, C(i=10Ci, and Cg=G.lC\. This region is shown 
by the shaded area which satisfies Ai5o,o,-e,o,i,o < 0, 

e,0,0,-e,0,0 < 0, ASq, -e, 0,0,0, 1 < 0, A_Eo, -e, -e, -e,l ,0 > 

0, and A£'o,o,o,+e,o,i > 0. The meaning of each line is 
as follows: 1: AiJo,o,-e,o.i,o=0, 2: Ai5_e,o,o,-e,o,o=0, 3: 

Ai5o,-e, 0,0, 0,1=0, 4: AiJo,-e,-e,-e,l,O=0, 5: A_Eo,0,0, + e, 0,1=0, 

where the expression is given by Eq. (|59p . 
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